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Abstract 

We numerically analyze recent high energy heavy-ion collision experiments based 
on a hydrodynamical model with phase transition and discuss a systematic change of 
initial state of QGP-fluid depending on colliding-nuclei’s mass. 

In a previous paper, we formulated a (3+l)-dimensional hydrodynamical model for 
quark-gluon plasma with phase transition and discussed numerically the space-time 
evolution in detail. We here compare the numerical solution with the hadronic distri¬ 
butions given by CERN WA80 and NA35. Systematic analyses of the experiments with 
various colliding nuclei enable us to discuss the dependences of the initial parameters 
of the hydrodynamical model on colliding nuclei’s mass. Furthermore, extrapolating 
the present experiments, we derive the possible hadronic distributions for lead-lead 
150GeV/A collision. 
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1 Introduction 


Physics of quark-gluon plasma is one of the most important problems in high energy physics. 
Quark-gluon plasma state is expected to be produced in extremely high temperature, which 
will be realized in ultra-relativistic heavy ion collisions. To obtain the higher temperature, 
experimental setting has been growing in three directions: the higher incident energy par 
nucleon, the heavier incident nucleus, and the heavier target nucleus. In any case, colliding 
energy is getting larger, but the larger energy does not directly mean the higher temperature. 
It is not obvious, whether large fraction of colliding energy will be transferred into thermal 
energy (temperature) or kinetic energy (systematic flow). Especially, the heavier target does 
not correspond to the larger colliding velocity in the center of mass system. Hence, if we 
stand on the Landau-type “simple stopping” picture, achieved energy density must become 
lower for larger target mass. In this paper, we will analyze the recent experimental data for 
various target nuclei and discuss systematic change of initial parameters in hydrodynamical 
model. 

In a previous paper [|l| we formulated a semi-phenomenological quantum transport theory 
for a quark-gluon plasma fluid based on an operator valued Langevin equation. Putting a 
mode spectrum and a damping into this formula, we can easily evaluate thermodynamical 
quantities and transport coefficients. Introducing a simple model spectrum to this formula, 
we have already discussed the space-time evolution of (1+1)-dimensional viscous quark- 
gluon plasma fluid 0 and (3+l)-dimensional perfect fluid quark-gluon plasma with phase- 
transition [Q. We have also discussed the case of baryon-rich quark-gluon plasma with phase- 
transition [|4j, in which one of the conclusions was that the baryon-number effect was not 
so much large. Hence, although the existing experiments do not seem to be baryon-free, we 
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apply a simple baryon-free quark-gluon plasma model to analyses of available experiments. 


2 Hydrodynamical model with phase-transition 


Following P|, we here introduce a simple model mode spectrum, 



1 + tanh T , Tc 

k-— 


( 1 ) 


where A is a parameter that adjusts the degrees of modes and Tq is the critical temperature 
responsible for the reliant phase transition. Supposing that the high temperature phase is 


dominated by massless u-, d- , s-quarks and gluons and the hadronic phase is dominated by 
pions and kaons, we can put A = 1.89, M = 200 MeV, d — 2 MeV, and Tq = 160 MeV. 
With these parameters, we can easily obtain an equation of state with phase transition-like 
behavior (fig. 1), which seems to reproduce the Lattice QCD result 0. 

The hydrodynamical equation is given by, 


( 2 ) 


d tl T ,w = 0, 


= EU fl U u — P(g^ v — U^U U ), (3) 

for perfect fluid. Here E,P and £/ M are, respectively, energy density, pressure and local four 
velocity. For cyliudrically symmetric expansion along the collision axis, it is convenient to 
introduce new variables, r, ?/, r and (p defined by 


t = t cosh ?/, 


r sinh r/, 


x 


r cos 0, 


y = r sm 0, 


3 







instead of ordinary coordinate t, x, y and z. As for the local four velocity, we represent it 
by four components, U T , U v , U r and associated with the new variables, 

U l = U T cosh rj + U v sinli 77 , 

U z = U T sinli q + U v cosh r/, 

U x = U r cos cf> + sin </>, 

U y = U r sin 0 — coscj). 

Because of = 1 and the cylindrical symmetry, the four components can be reduced to 

two variables Yt and Yl 0 , 

U T = cosh F t cosh(F L — r/), 

U 71 = coshF T sinh(F L — rj), 

U r = sinhF T , 

U* = 0. 

Taking account of the formation time, let us put the temperature distribution on rj and 
r at r = r 0 , 

(Ini — /7n ) 2 (t — r n ) 2 

T(r 0 , V ) = To^M- ^! 2 V a y (\rj\ - Vo ) - ~ r 0 )), (4) 

assuming that the hydrodynamical expansion starts at r = To = 1 fm later than the collision 
instance at r = 0 fm. To represents the initial temperature, y o and ro are measures of the 
longitudinal and transverse spreads in 77 and r, respectively. T 0 , rj 0 , r 0} a v and o> are input 
parameters to characterize our model on which we should impose a constraint given by a 
fixed value of the initial fluid energy E ini = x^tot ( E tot standing for the total collision energy 
and x f° r th e inelasticity). As for the initial condition of the local velocity, we use Bjorken’s 
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scaling solution, Yl = r), in the longitudinal direction, by which we can take the initial 
longitudinal flow into account. Initial values of these parameters should be determined on 
the basis of an appropriate physical discussion, but details of these parameters are deeply 
connected to the interaction mechanism and the thermalization process. In this paper, 
standing on phenomenological point of view, we try to choose these parameters so as to 
reproduce experimental results of both momentum spectrums. 

3 Particle Distributions 


The numerical solution of the hydrodynamical equation gives us the momentum distribution 
of hadrons, coming out from a local system with volume da M in local equilibrium with freeze- 
out temperature Tf, through the formula 

d 3 AN 1 U u da u ^ 

dp 3 VP 2 + m ‘ 2 exp - 1 (2vr) 3 

where p° = \Jp 2 + m 2 with hadron mass m (see ref. [Q). Integrating eq.(^) on the hypersur¬ 
face with Tf, we obtain 

dN r U^p^ U u da u P-£dP^dp , . 

d¥ = J exp(^p)-l (X ) 3 U 

for the Y-distribution, 


dN j' 1 1 U u da u dP^dp 

dr]' J VP 2 + m 2 exp (^p-) - 1 (2vr) 3 

for //-distribution, where rf = |ln stands for the pseudorapidity distribution, and 

dN r IT'Pn 1 U u da u dPi J dp 

PadPa = J VP 2 + ™ 2 exp(^)-l M 3 

for the Pr-distribut ion. 
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Our numerical results given by eq.(Q) is shown in fig.2. As is well known, the freeze-out 
temperature Tj is directly connected with the P? slope and we should determine it from the 
experimental Pj-spectrum. However, throughout this paper, we fix it equal to Tj = 140MeV 
because of the lack of sufficient experimental data. 


Our numerical solution of hydrodynamical model is so designed as to describe the head- 
on collision of the same kind of nuclei in the center of mass system. In order to apply 
the solution to the asymmetric collision, such as WA80 and NA35, we estimated effective 
participants in large target nucleus as, 


A eff = A- 


tar 


1 - 



A 


proj 


A 


tar 


3-i 

2 


( 9 ) 


A tar and A pro j being target mass and projectile mass, respectively. 


4 Numerical results 


Our hydrodynamical model contains many parameters to be adjusted. We have already 
fixed freeze-out temperature as Tf — 140MeV from fig.2 in the previous section. We may 
put formation time equal to Tq = 1 fm and initial transverse size equal to the size of smaller 
nucleus, 


r 0 = 1.2 x tfA proj - a r (fm), (10) 

where a r — 1 (fm). As for the local velocity, we use Bjorken’s scaling solution, Y L = rj, 

neglecting initial transverse flow. Therefore, the residual parameters to be adjusted are only 
initial temperature To at central point , and longitudinal size at initial time, r/ 0 and a v . 
We try to adjust these parameters so as to reproduce existing experiments. Our results for 
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sulphur beam experiments are shown in fig.3, and the values of parameters for sulphur and 
for oxygen beam are summarized in table 1. In any case, our model seems to reproduce 
the experimental results well with plausible values of the parameters. Table 2 stands for 
the results given by a equation of state in the limit Tq —> oo, which corresponds to a hot 
hadron-gas model. In this case, x became an unphysical value, larger than unity. It means 
that we failed to reproduce experimental results with the hot hadron-gas model. 

These results tell us that larger nuclei are more effective to make a high temperature 
fluid than smaller nuclei can do even in target experiment. For larger targets, the initial 
temperature must become higher, the longitudinal size smaller, and the inelasticity larger. 
This kind of tendency is just desired for the experimental quark-gluon plasma production. 

Sometimes, it has been pointed out that the present CERN energy is not high enough 
to make the Bjorken’s scaling quark-gluon plasma, and the stopping picture of Landau type 
is then realized. However, the simple stopping picture is not in agreement with our results. 
In the simple stopping picture, the only collision parameter is the colliding relative velocity 
7 , and the longitudinal size contracts with l/q and the energy density increases with q 2 . I 11 
target experiments, the larger target mass means the higher collision energy, but does not 
means the higher collision velocity. Table 3 shows the our results and relative velocities. The 
tendency is opposite to the simple Landau picture. Of course, it only means the failure of 
the simple stopping picture and we may describe the phenomena which require to introduce 
additional phenomenological parameters, such as energy dependent inelasticity 0. 

In this paper, we have fitted hydrodynamical parameters in a purely phenomenological 
manner without resorting to any dynamical model for initial collision process. Although we 
need a dynamical model for the quark-gluon plasma formation to understand the results in 
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the table 1, we may extrapolate our phenomenological results in order to obtain the rough 
sketch of coming experiments. Figure 4 indicates us that parameters in Pb + Pb 150 GeV/A 
collision are T 0 = 250MeV and rj 0 + a v = 0.8. The possible psuedorapidity-distribution of 
hadrons is shown in fig. 5 and the total charged hadronic multiplicity (pions and kaons) 
reaches 2095. 


5 Concluding Remarks 

We have applied our (3+l)-dimensional hydrodynamical model with phase-transition to the 
recent heavy-ion experiments. Our results show that the inelasticity increases with target 
mass, the longitudinal size is getting smaller with target mass, and the initial temperature 
becomes higher. Such a tendency is very promising for our aim of producing quark-gluon 
plasma experimentally and also for understanding the future experimental results at CERN 
and RHIC. 

The authors are indebted to Professor I. Ohba for his suggestive discussions. 
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Table 1. Phase-Transition Model 


Collision 

A 


Vo 

To(MeV) 

X 

S + A1 b 

27 f 

1.5 

1.0 

169 

0.52 

s + s 

32 

1.5 

1.0 

167 

0.46 

S + Cu ^ 

64 

1.5 

0.8 

180 

0.69 

S + Ag 0 

107 

1.5 

0.7 

188 

0.75 

S + w 

184 

1.2 

0.5 

202 

0.68 

S + An h 

197 

1.5 

0.7 

195 

0.88 

0 + Cu 

64 

167 

3.0 

0.3 

0.38 

0 + Ag 4 1 

107 

172 

3.0 

0.2 

0.54 

0 + An 4 ) 

197 

176 

3.0 

0.1 

0.62 

Data 1) from ||, 2) from 

Table 2. H 

|, 3) from [^] and 4) from |J 

adron Model 

Target 

A 

dr/j 

Vo 

T 0 (MeV) 

X 

Al 4) 

27 f 

0.8 

0.5 

335 

0.82 

S 2 > 

32 

0.6 

1.0 

312 

0.93 

Cuh 

64 

1.0 

0.5 

363 

1.37 

Ag lj 

107 

1.0 

0.5 

380 

1.49 

W 

184 

0.8 

0.7 

405 

1.65 

Au 4) 

197 

0.9 

0.7 

400 

1.77 


Data 1) from ||, 2) from |f and 3) from [0 


Table 3. Parameters for Landau Picture 


Collision 

Effective 

Participants 

Elab 

(GeV) 

Ecm 

(GeV) 

'7cm 

v v + Vo 

S + A1 

30.9 + 27 

6175.72 

561.218 

11 

1.5 + 1.0 

S + Cu 

32 + 49.6 

6400 

774.461 

8.3 

1.5 + 0.8 

S + Ag 

32 + 63.0 

6400 

873.332 

7.4 

1.5 + 0.7 

S + An 

32 + 81.1 

6400 

990.954 

6.5 

1.5 + 0.7 

O + Cu 

16 + 34.0 

3200 

453.783 

7.1 

3.0 + 0.3 

O + Ag 

16 + 41.9 

3200 

503.616 

6.4 

3.0 + 0.2 

O + Au 

16 + 52.7 

3200 

565.62 

5.7 

3.0 + 0.1 
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Figure Caption 


Fig-1 

The energy density as a f un ction of temperature; the solid line stands for our phase-transition 
model, the dotted line for the simple hadron-gas of pions and kaons, and the dot-dashed line 
for the quark-gluon plasma of massless u-, d-, s-quarks and gluons. In the phase-transition 
model, phase-transition takes place at Tc = 160MeV in the width of d = 2MeV. 

Fig.2 

The Pt spectrum of pions in S+S 200 GeV/A collision together with the data of NA35. The 
solid line stands for Tf = 140 MeV, the dashed line for Tj = 120MeV, the dotted line for 
Tf = lOOMeV, and dot-dashed line for Tf — 80MeV, respectively. 

Fig.3 

The pseudorapidity distributions of charged hadron in 200GeV/A Sulphur beam collision. 
The solid line stands for S+Al, the dashed line for S+Cu, the dotted line for S+Ag, and the 
dot-dashed line stands for S+Au, respectively. In our calculation, we take account of pions 
and kaons only. Plots are the data from ||. In the case of S+Al, ro is determined from the 
size of Al. 

Fig.4 

The values of initial parameters as a function of Total Collision Energy. The solid line 
stands for initial temperature,T 0 , and the dashed line for longitudinal extent, rj 0 + a v . The 
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plots from the left to the right correspond to O+Cu, O+Ag, O+Au, S+Al, S+Cu, S+Ag, 
and S+An, respectively. Plots at the right end are the extrapolation of parameters to Pb+Pb 
150 GeV/A collision. 

Fig.5 

The psuedorapidity distribution of charged hadrons (7r ± and K ± ) in Pb+Pb 150 GeV/A 
collision. Here, we put T 0 =250MeV, rjo =0.2, and a v = 0.6. 
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